Selection of den sites and chronology of denning by black bears in the eastern Sierra Nevada and western Great Basin

Abstract Use of dens during winter is an important strategy for American black bears (Ursus americanus) for both energy conservation and reproduction; and occupancy of suitable den sites has implications for reproductive fitness. Denning strategies may change as a result of changing climatic conditions and habitat loss. Black bears occupy arid environments in the eastern Sierra Nevada and the western ranges of the Great Basin Ecosystem. Our objectives were to identify: (1) which physical characteristics of habitat influenced selection of den sites at multiple spatial scales and (2) which environmental factors influenced timing of entrance and exit of dens by females and males. We evaluated selection of den sites by black bears at three spatial scales (300, 1000, and 4000 m) from 2011 to 2022. Terrain ruggedness was important for selection of den sites at all spatial scales. Within a 300‐m buffer from the den, bears selected den sites with rugged terrain, lower horizontal visibility, and greater canopy cover, resulting in more concealment and protection than that of the surrounding environment. Within 1000‐ and 4000‐m buffers around each den, bears selected den sites with rugged terrain, northern aspects, and steep slopes. At the 4000‐m scale, we observed interactions between sex with slope and distance to roads; females selected den sites on steeper slopes and closer to roads than did males. Females remained in the dens longer than males by entering earlier in the autumn and exiting later in the spring. Male bears exited their dens earlier with increasing consecutive days above freezing temperatures, but that relationship was weak for females. Knowing what characteristics are important for selection of den sites, and influence timing of denning, will be important for understanding how shifting climatic patterns will affect bears, particularly in arid environments that may be prone to wider fluctuations in climatic drivers of denning in the future.


| INTRODUC TI ON
Multiple species that experience environments with seasonal fluctuations in weather and availability of food have evolved strategies to survive in those conditions.Those strategies include moving to more mild climate conditions though migration (Berger, 2004;Hayes, 1995), seasonal breeding (Fuglei & Ims, 2008), or through periods of inactivity using torpor or hibernation (Bieber et al., 2012;Fowler et al., 2021;Hellgren, 1998;Humphries et al., 2003).
Hibernation is described as suppression of body temperature and metabolism for extended periods of time, which is coincident with use of dens by many species including members of the family Ursidae (Fowler et al., 2021;Hellgren, 1998;Melvin & Andrews, 2009).Use of dens over winter is a strategy used by many species of mammals to survive seasons of extreme weather and food scarcity (Humphries et al., 2003).Broadly defined, a den is a secure area or hideout where an animal spends time and is primarily used for rest, sleep (including torpor or hibernation), or reproduction (Davis et al., 2012;Libel et al., 2011;Robitaille et al., 2020).Dens used by bears (Ursus spp.) are areas where bears spend the winter, and can range from a hollow tree, under a rock or pile of rocks, under a shrub, to a shallow depression in the ground, bears may excavate or modify their dens or use the existing structure (Bard & Cain III, 2020;Beecham et al., 1983;Pelton et al., 1980;Ryan & Vaughan, 2004).A den site describes characteristics in immediate area around the location of the den.
Denning is an energetically advantageous strategy when the net costs of building and using a den are lower than the cost of remaining active during the same time period (Fowler et al., 2021;Humphries et al., 2003).Dens are used to protect both adults and young; and individuals select and use dens for thermal stability, energy efficiency, shelter from weather, protection during parturition, and protection from inter-and intraspecific predation on adults and young.Dens may be located in close proximity to sources of food that would be available upon exit (Boutros et al., 2007;Laurenson, 1994).
All species of bears in North America, including American black bears (Ursus americanus), polar bears (Ursus maritimus), and brown bears (Ursus arctos) use dens during winter for either hibernation, parturition, or both (Hellgren, 1998;Nelson et al., 1980).The time period when bears are within dens is not only important for energy conservation, but also for reproduction because parturition occurs and females begin to care for offspring while occupying dens (Anderson et al., 2012;Bard & Cain III, 2020;Pelton et al., 1980).
Bears cannot readily flee dens if disturbed without substantial energetic costs, and disturbance also can result in lowered body condition of adults and reduced survival of young (Baldwin & Bender, 2008).
Dens that were more energetically efficient resulted in brown bears emerging from dens in better body condition than those occupying dens that were less energy efficient (Shiratsuru et al., 2020).Time spent in dens during winter may vary across a species range or annually with fluctuations in climate and availability of resources (Amspacher et al., 2023).Therefore, selection of an appropriate den site is important for population growth or stability as well as survival and recruitment of young.Previous work has suggested that the type of den that bears selected was consistent with the type of cover that was most common in the surrounding area (Baldwin & Bender, 2008;Bard & Cain III, 2020;Beecham et al., 1983;Fowler et al., 2019;Schafer et al., 2018).
Timing of entry and exit from dens is also an important aspect of den occupation, and could affect both body condition of bears and reproductive success (López-Alfaro et al., 2013).Previous work has demonstrated differences in chronology of den entrance and exit between females and males; generally, males enter dens later and remain in the den for a shorter amount of time than do females (Beckmann & Berger, 2003;Fowler et al., 2019;Waller et al., 2012).Females have higher reproductive demand than males because of pregnancy and parturition, and therefore may enter the den earlier and exit later because of those demands (Doan-Crider & Hellgren, 1996;Friebe et al., 2014).Johnson et al. (2017) reported that female bears with dependent young entered dens earlier than females without young.Body condition of female bears can influence survival of young, so hibernation to conserve energy until parturition is especially important for females (Noyce & Garshelis, 1994).
Previous work has found timing of entrance and exit of dens to be variable (Bard & Cain III, 2020;Beecham et al., 1983;Fowler et al., 2019;Waller et al., 2012), and the period of hibernation can range from 0 to 212 days depending on seasonal food availability and winter severity (Fowler et al., 2021).Timing of entrance specifically has been correlated with decreasing abundance of forage plants (Fowler et al., 2019;Pigeon, Stenhouse, & Cote, 2016).While a specific environmental trigger or initiation of den entrance is not fully known, lower autumn temperatures, lower food availability, presence of snow, and incoming storms have been correlated with bears entering dens (Baldwin & Bender, 2010;Fowler et al., 2019;Friebe et al., 2001;Schooley et al., 1994).Food subsidies from areas of human habitation may also affect denning by providing consistent food sources throughout the year, even causing some individuals to forego denning entirely (Amspacher et al., 2023;Krofel et al., 2017).
Previous research has also shown temperature to influence exit from dens, and movement away from dens following emergence may in part be driven by warming temperatures (Fowler et al., 2019;Gonzalez-Bernardo, Giulia, et al., 2020;Johnson et al., 2017;Waller et al., 2012).
The climate, aridity, and topography of the eastern Sierra Nevada and western Great Basin as well as proximity to humans may cause patterns of denning to differ from other areas where denning behavior had previously been studied.The objective of this study was to investigate selection of den sites and timing of entrance into and exit from dens by black bears in the eastern Sierra Nevada and western Great Basin of Nevada.We hypothesized that bears would select primarily rock or tree dens across our study areas, because large trees are common in the Sierra Nevada and rock piles are common in both the Sierra Nevada and Western Great Basin.We hypothesized that den sites would have high concealment and would be located far from human impacted areas and sources of disturbance.
Based on that hypothesis, we predicted that bears would select dens with higher concealment (low visibility) and greater distance to the nearest road than was generally available in the study area.
Additionally, we hypothesized that females and males would differ in type of den selected, such that females would select den types that provide more protection than dens selected by males.In investigating chronology, we hypothesized that mean entrance and exit dates would differ among sexes and that females would enter dens earlier and exit later than males.We hypothesized that entrance and exit would be correlated with temperature and the presence or absence of snow preceding entrance or exit from dens.We predicted that earlier entrance dates within the dataset would be correlated with cooler temperatures and deeper snow in autumn, and that earlier exit dates would be correlated with warmer spring temperatures and lower snow depth.

| Study area
The study area consists of two eco-regions in northwest Nevada, USA: the Sierra Nevada Mountains east of Lake Tahoe (eastern Sierra Nevada) and the western Great Basin (Figure 1).This study area is unique in that it is one of the driest regions occupied by black bears, and there is a large wildland-urban interface, meaning bears may forage in human occupied areas, but den in wildlands (Beckmann & Berger, 2003;van Manen et al., 2019).While black bears may now be found across much of Nevada (Lackey et al., 2013), our study area includes the mountain regions where bears are most common: The Carson Range in the eastern Sierra Nevada, the Pine Nut Mountains, and the Virginia Range in the western Great Basin, as well as surrounding areas and nearby basins (Figure 1).Elevations range from a maximum of 3316 m in the Carson Range, 2397 m in the Virginia Range, and 2882 m in the Pine Nut Mountains, to minimum elevations within the basins around 1200 m.

| Data collection
Our dataset consisted of den-site locations of black bears from 2011 to 2022, which was provided by the Nevada Department of Wildlife (NDOW).A total of 116 den sites were identified over that time frame and NDOW biologists visited 76 of those sites from 2011 to 2022 to confirm location, presence and sex of the bear, and type of den.Locations of den sites were determined by NDOW using GPS (Global Positioning System) collared bears and were identified as a den either when GPS locations formed a cluster of spatial points over the winter, when satellite transmissions or loss of reception from the collar indicated that a bear had quit moving, or data from spring indicated that the bear had exited the den.We removed five dens from our study because they occurred underneath human homes, and bears were encouraged to vacate the structure by NDOW biologists.
We grouped dens into five types: "rock," "tree," "excavated," "exposed," or "other."Bears in our study either used an existing structure such as a cavity under rocks or inside a hollow tree, but there were also indications that bears modified those sites or excavated new sites.We defined rock dens as a den underneath a single rock or pile of rocks with a cavity below.We defined tree dens as a standing tree that was dug out underneath by the bear or was already hollow.Dens excavated by bears had no preexisting structure, but were made up of shallow holes in the ground or used brush piles.Exposed dens were areas with no excavation around the den site and no direct use of surrounding trees or rocks, although exposed dens were often found at the base of trees.We also included a group of five dens, which were categorized as "other."Those dens occurred in anthropogenic structures (not houses), which included culverts under roads and abandoned mines.
In addition to the 76 dens visited by NDOW biologists from 2011 to 2022, we visited 26 additional den sites to collect data on vegetation and habitat structure, predominantly in the Carson Range of the Sierra Nevada from August to December 2022.At those 26 sites, we measured characteristics from the den opening within a 15-m diameter circular plot around the den and two randomly generated sites within a 300-m radius from the den (Bard & Cain III, 2020).
Those characteristics included: dominant vegetation type, ground cover (%), number of trees, tree cover (%), horizontal visibility (%), and diameter at breast height (DBH) of the two largest trees to generally characterize the size of trees at each site (Tables 1 and 2).The 26 dens selected for field data collection were randomly selected from the full set of 116 dens, but we also considered feasibility and safety of accessing those sites.We visually estimated the dominant vegetation type as one of the following categories: bare ground, shrub (denoting family of shrub), trees, or leaf litter.We also visually estimated percent cover of those vegetation categories, including shrub species.For measurements involving estimations, the same researcher made the estimates each time to ensure consistency of measurements among sites.We measured tree cover (%) at five sites within each plot using a spherical convex densiometer (Lemon, 1956) at the center of the plot facing north and then at the edge of the plot along each of the four cardinal directions facing the center (Bard & Cain III, 2020;Pigeon, Cote, & Stenhouse, 2016).We measured horizontal visibility using a 1-m tall cylinder with a 33-cm radius (Bard & Cain III, 2020;Ordiz et al., 2009).We placed the cylinder at the center of the plot and recorded the percentage of the cylinder that was visible from each cardinal direction at the edge of the plot observed from 1 m above the ground.We calculated a minimum distance of total concealment, by measuring the minimum distance away from the den in which the cylinder was no longer visible (Bard & Cain III, 2020).
In addition to measuring characteristics in the field at 26 den sites, we measured geospatial variables at the full set of 116 dens.
Each of the 116 den sites was paired with randomly located points that were generated using ArcGIS (Esri, Redlands, CA, USA).Those random points were generated within three different buffer zones around the den: a 300-, a 1000-, and a 4000-m radius.This paired design was modified from Pigeon, Cote, and Stenhouse (2016) who used random points within 1500 m of each den site and Bard and Cain III (2020) who used random points within 50-250 m from each den site.We included a third spatial scale at 4000 m based on previous work that has shown habitat selection for both black and grizzly bears to be influenced by features up to 4000 m away (Beckmann et al., 2015;Mattson et al., 1986).Using those three different buffer zones allowed us to investigate selection at both local and broad spatial scales (Gray et al., 2017).
To ensure that we were appropriately characterizing availability from each of our buffer zones, we tested our models with iteratively increasing numbers of random points paired with each den site.For TA B L E 1 Summary statistics, means ± standard deviations, of geospatial landscape-level characteristics and local characteristics measured in the field at den sites of black bears (n = 26) and random sites within a 300-m buffer (n = 52), and remotely at den sites (n = 116) and random sites (n = 9280, 17,400) within 300-, 1000and 4000-m buffers.We extracted geospatial variables including values of slope, tree cover, elevation, distance to the nearest road, aspect, and terrain ruggedness at each den site and random location (Table S1, Table 1).
Slope, elevation, and aspect were extracted from a digital elevation model (30-m resolution).Tree cover (30-m resolution) and distance from each den and random location to the nearest road (30-m resolution) were both acquired from the Western USA Geodatabase at the University of Nevada Reno.Our estimate of ruggedness is a modified version of the vector ruggedness model (30-m resolution) (Dilts et al., 2023;Sappington et al., 2007).
To investigate entrance and exit dates based on environmental covariates, we obtained environmental data from Snowpack Telemetry (SNOTEL) sites managed by the USDA Natural Resources Conservation Service (NCRS) National Water and Climate Center.
We chose four SNOTEL automated data collection sites, Mt. Rose (Site 652), Heavenly (Site 518), Marlette Lake (Site 615), and Tahoe City (Site 809), as being representative of the den-site locations in our study (mean distance between the den and the nearest SNOTEL site was 9.2 ± 6.49 km); and for each environmental covariate, we used mean values of the four sites at each day (Figure 1).Entrance date was determined by NDOW biologists as the date where collar locations become stationary in the winter, or the last transmission of that collar for the winter.Exit date was determined by NDOW biologists as the date that the bear left the den site and continued traveling away from the site for two consecutive GPS locations without returning.Some den sites did not have entrance and exit dates recorded, so we used a subset of dens for each analysis: entrance date (n = 68) and exit date (n = 42).We also excluded dens that were too far from a SNOTEL location to document snow depth or temperature.For example, we excluded den sites from the Pine Nut Mountains and Virginia Range in the western Great Basin because there are no SNOTEL sites in those ranges.Our chosen environmental covariates were minimum daily temperature and daily snow depth.We also included elevation in the entrance and exit models because of differences in amounts of snow accumulation by elevation, which was extracted from ArcGIS.In order to test a window of time that would be ecologically meaningful based on the likely time frame that bears were preparing to enter dens, we averaged the daily snow depths of the 2 weeks preceding entrance or exit for each den site (Gonzalez-Bernardo, Russo, et al., 2020).To determine the effects of temperature on entrance and exit from dens, we used the number of consecutive nights below freezing for entrance and above freezing for exit.

| Statistical analysis
We began all analyses by standardizing variables using a ztransformation (Zar, 2010).We transformed aspect, which is a circular variable, using both a sine (east-west) and a cosine (northsouth) transformation, prior to standardizing (McKee et al., 2015;Zar, 2010).We then assessed collinearity for all predictor variables using a Pearson correlation matrix in R (4.2.2 R Core Team, 2022) and based on biological relevance, eliminated one of any two variables in the same model that were highly correlated |r| > .70 with one another (Heffelfinger et al., 2020;Long et al., 2014;Stewart et al., 2002).In our analysis, no variables were highly correlated.
Since we were interested in investigating multiple scales of selection, we analyzed each of the four datasets separately, resulting in To quantify characteristics of den sites selected by bears, we used conditional logistic regression using the "clogit" function in the survival package in R version 4.2.2 (R Core Team, 2022; Therneau, 2022;Therneau & Grambsch, 2000).Conditional logistic regression allowed us to use a matched-case control format for used and available (random) locations (Bard & Cain III, 2020;Manly et al., 2002;Schafer et al., 2018).We used the "stepAIC" function in the MASS package in R (Venables & Ripley, 2002) to rank models by Akaike's Information Criterion corrected for small sample size (AICc; Burnham & Anderson, 2002) and to identify a top-performing model.We then evaluated our top-performing model for uninformative effect sizes defined as those with 95% intervals that crossed zero (Arnold, 2010).Because of small sample sizes, we were unable to include den type in our analyses, but we included interactions with covariates with sex to address differences in selection of den sites between males and females.
We assessed differences in types of dens between females and males using a Fisher's Exact Test (Zar, 2010).We then used a series of two-sample z-tests for proportions on each type of den to look for differences between males and females within each den type.
For the z-tests, we adjusted α from 0.05 to 0.01 by dividing 0.05 by the number of z-tests performed (5) using a Bonferroni correction for repeated comparisons (Zar, 2010).We tested for regional differences in den locations and types between the Sierra Nevada (Carson Range) and the Great Basin (the Virginia Range and the Pine Nut Mountains) using Fisher's Exact Test and a series of two-sample z-tests.
We used a two-sample t-test to test for differences in den entry and exit dates based on sex, after converting entry and exit dates to ordinal dates and adjusted entry dates for those that overlapped January 1 (Zar, 2010).We used one-way analysis of variance (ANOVA) on means of den entry and exit dates to test for differences among years (Zar, 2010).In 2013 and 2018, there were no exit dates recorded, so we were unable to include those years in the analysis.
To test the effect of static and time-varying environmental conditions on individual dates of entrance and exit from den, we fit denning chronology models (analogous to known-fate survival models) in a Bayesian framework using JAGS (Just Another Gibbs Sampler) (Plummer, 2003).In the "entrance" model, bears entered den with a daily probability p ent , and in the "exit" model, bears emerged from den with a daily probability p exit .Both p ent and p exit were modeled as a logit-linear function of covariates including elevation, snowpack, consecutive days with sub-freezing temperatures, day of year, and sex (see below for details).Entry and exit from den in the exit and entrance models, respectively, was assumed to represent an absorbing state (e.g., once a bear had entered den in the "entrance" model, we assumed it could not exit).Because we were interested in time periods two weeks before beginning hibernation, we quantified the number of consecutive sub-freezing days (≤−10°C) prior to each potential entrance date ("consecutive freeze days" for modeling the timing of entrance to den) and the number of consecutive above-freezing days (>0°C) prior to each potential exit date ("consecutive warm days" for modeling the timing of exit from den).Similarly, we used the SNOTEL data to quantify the mean daily snow depth across the 2 weeks before each potential entrance and exit date (Table S1).We also included the (static) den elevation in these models because of variation in temperature and snow depth by elevation.Because we hypothesized differences in overall timing patterns between female and male bears, we included sex as a covariate in both the entrance and exit models (males were treated as the reference class).Finally, we tested for an interaction between sex and the number of consecutive freeze days (entrance model) or warm days (exit model).To visualize the effects of covariates (e.g., elevation, number of consecutive sub-freezing days) on denning chronology, we used our models to derive the median date of entry and exit (50% quantile of the distribution of entry or exit dates) across specific gradients (holding all other covariates at their mean values).
We assigned uninformative uniform priors (min = 0, max = 1) on all probabilities and weakly regularized Gaussian priors (mean = 0, var = 10) to each regression coefficient.We ran each model with three chains and 10,000 iterations, and assessed convergence using visual assessment of trace plots and ensuring all potential scale reduction factors (R-hat) were below 1.1 (Brooks & Gelman, 1998).

| RE SULTS
Our hypothesis that females and males differed overall in types of dens used was weakly supported (Fisher's exact test p = .045);however, no pairwise comparisons (z-proportion tests) were interpretable after applying a Bonferroni correction (revised α of 0.01) (Figure 2).We also investigated proportional differences in use of to use tree dens more in the Sierra Nevada (n = 13) than in the Great Basin (n = 0, p = .04)(Figure 2).  1, Figure 3).
Probability of selection of den sites at 300-m scale declined with increasing horizontal visibility above ~50% (Figure 3).Probability of selection of den site was greatest in rugged terrain with more than ~40% tree cover (Figure 3).Included in the study were 70 individuals that had recorded either exit or entrance date from the den.In general, females entered dens earlier (p = .04)and exited later (p = .002)than did males (Figure 5).The mean date of entry for all bears across all years was December 6th (±17) (x ± SD) for females and December 16th (±18) for males.The mean date of exit for all bears across all years was March 28th (±16) for females and March 12th (±13) for males.When testing for differences among years, there was no difference in entry date for females (p = .7)or males (p = .6)among years (2013-2022) (Figure 5).There was also no difference in exit date for females (p = .6)or males (p = .8)among years (2014-2022, excluding 2018) (Figure 5).
Our results indicated a weak influence of the number of consecutive sub-freezing days on the daily probability of den entrance (and consequently on the median date of den entrance; Table 2, Figure 6), and a strong influence of the number of consecutive warm (above-freezing) days on the daily probability of den exit (Table 2), with strong consequences for the median date of exit; (Figures 6 and   7).We detected no temperature-by-sex interaction on the probability of den entrance, but there was evidence for a temperatureby-sex interaction on the probability of den exit, with males more responsive to the number of consecutive above-freezing days than females (Figure 7, expected consequences for the median date of entry illustrated in this figure).Mean snowpack was not strongly predictive for either den entrance or exit (Table 2).Elevation was weakly predictive of both the probability of den entrance and the probability of den exit (Table 2; Figures 6 and 7).

| DISCUSS ION
Bears occupied dens predominantly under rocks and within hollow trees, reflecting selection for hidden and protected dens.A large rock or tree would obscure much or all of the hibernating bear within the den cavity, suggesting that protection from disturbance was an important factor in choosing a den site.Although, Waller et al. ( 2012) observed bears using mostly tree and ground dens, bears in our study occupied some ground and exposed dens, but rock piles in both study areas and hollow trees in the Sierra Nevada were much more commonly used by bears in our study.An enclosed den may also offer environmental protection with stability in temperature and humidity that would not be present in an exposed or ground den.Pinyon pine and juniper trees, characteristic of the Western Great Basin, are smaller and less likely to be available as hollow trees for bears to use as dens compared to the large species of pines that are common in the Carson Range.We observed no tree dens used by bears in the Great Basin, and we were unable to document the presence of any hollow trees that could have been occupied as dens by black bears.Black bears inhabiting that region used primarily rock dens, which were readily available in that area.Goodrich and Berger (1994) reported a similar difference, recording 16 tree dens in the Carson Range in the Sierra Nevada, but only two in the Sweetwater Range in the Great Basin.
Spatial scale strongly influenced selection of den sites, and our model results varied with spatial scale, with the exception that ruggedness of terrain entered every one of our models.At all three spatial scales, including the 300-m model that consisted of both field collected and geospatial data, we observed dens to be located in more rugged terrain than was generally available.At the 1000-and 4000-m scales, den sites were located on steeper slopes and more northerly aspects than was available.Although the ruggedness metric is calculated from both slope and aspect, it is not correlated with those covariates (Dilts et al., 2023;Sappington et al., 2007).In the 300-m scale analysis that included field collected and remotely sensed data, characteristics that had the most influence on selection of den sites were canopy cover of trees, low horizontal visibility (e.g., higher concealment), and ruggedness of terrain, which supported our hypothesis that bears selected for high concealment of Only in the 4000-m buffered area were interactions of sex with slope and distance to roads detected.
through survival of offspring and body condition of individuals over winter (Baldwin & Bender, 2008;Shiratsuru et al., 2020).Our data indicated that including field-based measurements with the remotely sensed data at the smallest scale resulted in a model that showed strong selection for concealment of the den site in addition to being located in rugged terrain, which was the only variable that entered the remotely sensed model at that scale.Therefore, we suggest that at localized spatial scales including data collected in the field combined with remotely sensed data resulted in more effective understanding of selection of den sites.As we expanded the spatial scale that we modeled selection, the models became somewhat more complex.Ruggedness of the terrain in the area of a den site is likely the most important characteristic used by bears to select a den site since that characteristic entered both the field collected and all of the remotely sensed models.Ruggedness may add protection and concealment of the den site to reduce disturbance for hibernating bears (Hellgren & Vaughan, 1989;Hightower et al., 2002;Waller et al., 2012).Disturbance of bears in winter dens has been identified as a concern for den abandonment in our study system (Goodrich & Berger, 1994).
At the broad 4000-m scale, we detected interactions between sex with slope and sex with distance to roads.Females selected den sites on steeper slopes and closer to roads than did males.Excavated dens on slopes under 25° have been suggested to decrease thermal effectiveness of the den, in addition to increasing the risk of snow falling into the den opening or flooding in the spring (Crupi et al., 2020;Servheen & Klaver, 1983).The majority of dens used in our study were in hollow trees or under rocks, which provided more protection and may be less dependent on the effects of slope compared with excavated dens.Distance to the nearest road could influence disturbance, since dens that are further away from roads would be less likely to be disturbed by motor vehicles; males selected den sites farther from roads possibly to avoid those disturbances (Gantchoff et al., 2019;Lustig et al., 2021).Gantchoff et al. (2019) also observed female black bears denning closer to roads than males and suggested that this behavior by females, especially those with young, may be a way to reduce the risk of infanticide by males.We were not able to differentiate nonpregnant, pregnant, or females | that were accompanied by dependent young, and although other studies have noted that avoidance of males was less likely to be observed in nonpregnant females (Gantchoff et al., 2019) the signal was strong in our data.
One of the challenges of this study was that we used random sites to determine availability rather than "potential den sites."Humans perceive the world differently than do animals (Nams et al., 2006), and attempting to define a den site would lead to a much larger source of bias than using randomly generated locations.Therefore, we were unable to include den types in our models, and we also were unable to quantify availability of dens across the landscape.Bears in our study areas used multiple types of dens, including trees, rocks, exposed, or ground dens, and as a result almost any location could be a potential den site.
We observed almost no fidelity to or reuse of dens and saw only one instance of a den being reused, which was by a different individual.Similarly, other studies have reported few instances of individuals reusing dens, and den reuse is usually made by different individuals (Klenzedorf et al., 2002;LeCount, 1983;Linnell et al., 2000).Klenzedorf et al. (2002) not only reported some fidelity to den type in West Virginia, but also reported only a few instances of individuals reusing dens.Fidelity to a den was either absent or not recorded in many black bear studies, although there is some evidence of fidelity to den sites in areas with low availability of suitable dens (Davis et al., 2012).There is some evidence of fidelity to the same general area in brown bear populations (Sorum et al., 2019).Fidelity to a den has been hypothesized as an adaptation to low availability of existing den sites (Alt & Gruttadauria, 1984;Davis et al., 2012;Johnson & Pelton, 1980), which is unlikely to be a constraint in our study areas because there is an abundance of large trees and rock piles (Grayson, 2011).The lack of fidelity to dens that we observed suggests that suitable den sites are not limiting for bears in our study areas.Waller et al. (2012) noted that reuse of dens in the Southeast is typically in tree dens (Crook & Chamberlain, 2010;Linnell et al., 2000).Tree dens may be strongly selected in areas prone to flooding (Schwartz et al., 1987), which is not a selective factor in our arid study areas.
Similar to other studies of black bears, females occupied dens longer than males; both entering earlier and exiting later (Fowler et al., 2019;Mitchell et al., 2005;Waller et al., 2012), although that pattern did not occur in New Mexico (Inman et al., 2007).Timing of exit from dens was strongly related to the cumulative number of days with minimum temperatures above freezing for males, but that relationship was weak for females.Miller et al. (2017) reported that spring temperatures were negatively correlated with exit data and that when spring temperatures were lower, exit from dens was delayed.Doan-Crider and Hellgren (1996) suggested that emergence from dens by females was probably more closely related to development and growth of young than to food availability or weather conditions, which may explain why weather was less likely to affect exit of dens by females.Others have reported that reproductive state of females (not pregnant, pregnant, with cubs, or with yearlings) affected timing of exit and entrance to dens (Gantchoff et al., 2019;Immell et al., 2013;Johnson & Pelton, 1980).An important caveat to our analysis of den entrance and exit dates is that we were unable to include reproductive status of females or age of individuals, which has been shown to have as much if not more influence than environmental characteristics (Johnson et al., 2017).Pregnant female bears have also been seen to enter dens earlier than nonpregnant females (Fowler et al., 2019;Garshelis et al., 2020), thus lack of those data may be confounding our results.Timing of entrance and exit also has Mean snowpack had a very small effect size in this model and is not visualized here.
been tied to body condition, and bears in good body condition may den earlier than those in poor condition (Schooley et al., 1994).
We were unable to adequately quantify body condition of bears at the timing of entrance, nor were we able to document food availability, or to tie environmental conditions to potential food availability, but that could be an area of further research with this population of bears.
Our results regarding selection of den sites reinforce previous work, which reported that habitat selection varied relative to local versus landscape-level characteristics (Crook & Chamberlain, 2010;Reynolds-Hogland et al., 2007)

F
I G U R E 1 Map of 116 black bear den sites in the eastern Sierra Nevada (light blue points) and western Great Basin (dark blue points) identified from 2011 to 2022, and selected Snowpack Telemetry (SNOTEL) sites within the eastern Sierra Nevada.| 5 of 15 LONG et al.
four separate models: (1) both landscape characteristics and characteristics measured in the field compared with two random points close to the den (300-m), (2) landscape characteristics collected from geospatial data compared to 80 random points close to the den (300-m), (3) landscape characteristics compared with 80 random points at a middle distance from the den (1000-m), and (4) landscape characteristics compared with 150 random points at a far distance from the den (4000-m).
their dens.Interestingly, for the models that only included remotely sensed data, each spatial scale resulted in the addition of one or more characteristics with the 300-m model having only ruggedness in the model.Waller et al. (2012) have noted that evaluation selection of den sites is a hierarchical process and should be examined at multiple spatial scales to improve understanding of how bears select den sites.Additionally, factors that have the greatest effect on fitness | 9 of 15 LONG et al. in terms of selection should occur at greater spatial scales and lesser fitness limiting factors should occur at smaller spatial scales (McLoughlin et al., 2002; Rettie & Messier, 2000; Waller et al., 2012).Nevertheless, characteristics immediately around the den, such as concealment, and the structure and type of den affect energy efficiency and concealment, which also have direct effects on fitness, F I G U R E 3 Effect plots illustrating the functional relationship between the relative probability of den-site selection (Y-axis) and the three covariates included in the best performing (based on AICc) conditional logistic regression model for selection of black bear den sites identified within the eastern Sierra Nevada and western Great Basin from 2011 to 2022 (n = 26).Effect sizes combining geospatial and field collected characteristics (a).Curves of the relative probability of selection for horizontal visibility (b), tree cover (c), and terrain ruggedness (d) with 95% confidence intervals.F I G U R E 4 Curves of the relative probability of selection using conditional logistic regression for selection of black bear den sites in the eastern Sierra Nevada (2011-2022) (n = 116) using landscape-level geospatial characteristics.Spatial scales are: 300-m buffer with 80 random points per den (a), 1000-m buffer with 80 random points per den (b-d), and 4000-m buffer with 150 random points per den (e-h).

F
Ordinal date of entry (a) or exit (b) for male and female black bears in the eastern Sierra Nevada and western Great Basin from 2011 to 2022 (n = 70). of entrance was adjusted for entrance dates after January 1. Mean date of entry was December 10th for all individuals across all years (2011-2022).Mean date of exit was March 20th for all individuals across all years.Two-sample t-tests support that females enter dens earlier (p = .04)and exit later (p = .002)than males.F I G U R E 6 Cumulative probability of entry into den for each ordinal date during autumn for black bear dens in the eastern Sierra Nevada (2011-2022) (n = 68) based on a time series survival model (a).Mean snow depth 2 weeks prior to entrance (b), expected date of entry based on cumulative number of consecutive days below freezing (−5°C) prior to entrance (c), and elevation of den site (d).

F I G U R E 7
Cumulative probability of exit for number of days and nights above freezing for each ordinal date during spring for black bear dens (n = 42) in the eastern Sierra Nevada (2011-2022) based on a time series survival model (a).Expected date of exit based on day of the year, cumulative numbers of days above freezing 10 days prior to exit (b), elevation of the den (c), and elevation of den site (d).
. Nevertheless, our work provides some unique findings not shown in previous work.Our examination of selection at multiple spatial scales illuminates the importance of rugged terrain in selection of den sites.Additionally, cumulative days above freezing strongly affected timing of exit of dens for males, although that effect was weaker for females.Nevertheless, there may still be a knowledge gap in linking den-site selection specifically to reproductive fitness.Additionally, changes in the length and severity of winter affect food availability, which if connected to timing of den entrance and exit, could change bear behavior.Human encroachment into habitat could further affect selection of den sites with the potential of restricting availability of or access to den sites in wildland areas.